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CN ■ 

The quantum dynamics of an ensemble of interacting electrons in an array of random scatterers is 
treated using a new numerical approach for the calculation of average values of quantum operators 
and time correlation functions in the Wigner representation. The Fourier transform of the product 
of matrix elements of the dynamic propagators obeys an integral Wigner-Liouville-type equation. 
C/3 ' Initial conditions for this equation are given by the Fourier transform of the Wiener path integral 

representation of the matrix elements of the propagators at the chosen initial times. This approach 
combines both molecular dynamics and Monte Carlo methods and computes numerical traces and 
spectra of the relevant dynamical quantities such as momentum-momentum correlation functions and 
spatial dispersions. Considering as an application a system with fixed scatterers, the results clearly 
demonstrate that the many-particle interaction between the electrons leads to an enhancement of 
the conductivity and spatial dispersion compared to the noninteracting case. 
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q ■ I. INTRODUCTION 

Noninteracting electrons in an array of fixed random scatterers are known to experience Anderson localization at 
temperature T = in one-dimensional systems. In two-dimensional systems single-parameter scaling theory predicts 
that all states are localized as well Q. In three-dimensional arrays of random scatterers localization appears at 
the band extremities, while the center may consist of delocalized states, indicating a disorder driven metal-insulator 
■ transition in three dimensions only. 

This picture has been challenged by recent experiments, which suggest a metal-insulator transition also in disordered 
two-dimensional electron systems Q . A theoretical explanation || indicated that the Coulomb interaction between the 
electrons plays a central role in this effect. Moreover, the experimental study of persistent currents in mesoscopic metal 
rings ^ yields currents that are two orders of magnitude larger than predicted by theories based on noninteracting 
electrons ||. These findings suggest that it is the many-particle interaction which leads to derealization tendencies. 
To the best of our knowledge this has first been suggested by M. Pollak and coworkers || 0. 

Evidence for this influence of the Coulomb interaction has been obtained by examining the problem of two interacting 
electrons in a one-dimensional disordered band ||. It is, however, not clear whether this schematic model is able 
to yield answers for more realistic systems described by an ensemble of many electrons. Therefore, attempts have 
been made to calculate directly the conductivity of a disordered system of many electrons. Numerical results for the 
ensemble averaged logarithm of the conductance have been obtained for small systems by applying a Hartree-Fock 
based diagonalization scheme for electrons without |j| [l0| and with spin JTlJ. An increase of the conductivity for 
certain regions in the disorder-interaction parameter space has been demonstrated in these publications. 

The purpose of this paper is two-fold: i) We present our approach, which does not rely on small system sizes. 
This can be applied to a wide variety of different physical systems, such as plasmas, liquids, and solids, ii) As an 
illustration we study the problem described above. In particular, in this paper we investigate the influence of the 
many-particle interaction on electronic transport. We consider as a case study a one-dimensional disordered array of 
scatterers interacting repulsively with the electron system. Without electron-electron interaction such a system shows 
Anderson localization. It is the purpose of this application to study the effect of the long-range electron-electron 
Coulomb interaction on the mobility of the electrons. 

Anderson localization at temperature T = relies on quantum coherence of electron trajectories and results from 
interference. The key parameter in the physics of electron interference phenomena is the dephasing time of electrons. 
At finite temperatures the electron coherence is destroyed by both the electron-electron and the phonon-electron scat- 
tering [|l2| [l3| [l4|, |l5). To study the influence of these effects on kinetic electron properties in a random environment 
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we have simulated the quantum dynamics in a one-dimensional canonical ensemble at finite temperature for both 
interacting and noninteracting electrons using a Quantum-Dynamics-Monte-Carlo scheme. The main quantities cal- 
culated in this paper are the temporal momentum-momentum correlation functions, their frequency-domain Fourier 
transforms and the time dependence of the spatial dispersions. We discovered that the results strongly depend on the 
electron-electron interaction, clearly demonstrating the delocalizing influence of the many-particle interaction even at 
finite temperatures. 

Our approach also treats the positions of the scattering centers as dynamical variables. We are therefore in the 
position to generate various initial conditions, including the scenario of "electron bubbles" as an interesting side-result. 

We start in Section II by presenting the Quantum-Dynamics-Monte-Carlo approach used in this work, which is 
based of the Wigner representation of quantum statistical mechanics. It treats both electrons and scatterers on an 
equal footing. Details of our model, relevant correlation functions and the system of reduced units are introduced in 
Section III. Numerical results for the momentum-momentum correlation function in the time and frequency domains 
and for the spatial dispersion are given in Section IV. Finally, Section V contains a discussion of our results. 



II. WIGNER REPRESENTATION OF TIME CORRELATION FUNCTIONS 



According to the Kubo formula the conductivity is the Fourier transform of the current-current correlation function. 
Generally, time correlation functions Cfa^) — \F(fi)A(t)^ for a pair of dynamic variables F(t) — /p(t)J and 

A[t) = (A(t)\ are among the most interesting quantities in the study of the dynamics of electrons in disordered 

systems of scatterers, such as transport coefficients, chemical reaction rates, equilibrium and transient spectroscopy, 
etc. Our starting point is the general operator expression for the canonical ensemble-averaged time correlation function 



C FA {t) = Z- 1 Tr{Pe i6t ^ h Ae- i ^ R \ 



(1) 



where H is the Hamiltonian of the system expressed as a sum of the kinetic energy operator, K, and the potential 
energy operator, U. Time is taken to be a complex quantity, t c = t — ihf3/2, where (3 — 1/fcgT is the inverse 
temperature with ks denoting the Boltzmann constant. The operators F and A are quantum operators of the 
dynamic quantities under consideration and Z = Tr is the partition function. 

The Wigner representation of the time correlation function in a ^-dimensional space can be written as: 



C FA (t) = (2TrhY 



dnid^ F(ni) A(pi) W(pi; ^ t; ifr/3), 



(2) 



where we introduce the short-hand notation for the phase space point, fii — (pj, qi), (i — 1, 2), and p and q comprise 
the momenta and coordinates, respectively, of all particles in the system. In the definition (|2|), W{n\] [i2] t; ih/3) is 
the spectral density expressed as 



W(/ii;Ai 2 ;t;i/i/?) = Z~ 



d£,id£ 2 e K e R x 



jHt'/h 



and A(n) denotes Weyl's symbol [[TtJ of the operator A 
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and similarly for the operator F. Hence the problem of the numerical calculation of the canonically averaged time 
correlation function is reduced to the computation of the spectral density. 

To obtain the integral equation for W let us introduce the a pair of dynamic p, g-trajectories 
{Qt(t;pi, qi, t), Pt(t;pi, qi, t)} and {q T (T;p2,q2, t), p T (r;p 2 , 92, t)} starting at r = t from the initial condition {qi,px} 
and {q2,P2} propagating in 'negative' and 'positive' time direction, respectively: 



with 



dp _ 1 
~d7 ~ 2" 
Pt(r = i;pi,<7i,i) 



F[q T (T)] ; 



Pi 



<h = Pt(t) 

dr 2m ' 
Qt(r = t;px,qi,t) = q u 



(5) 
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and 



with 



Pt(r = t;pi,qi,t) = p x 



dq 



Pt{t) 



2m 

q t {r = t;p 1 ,q 1 ,t) = q u 



(6) 



while w{a,q) = * | dg' <7(g - g') sin ) - FU/)VA(.s,. 



(7) 



(5(s) is the Dirac delta function and F (q) = —VU with U being the total potential. 
Then as has been proved in JlJ, the W obeys the following integral equation: 

W(p, 1 ;p, 2 ;t;ih/3) = W(p , q ;Po, <?o; i^P) + 



~ J dr J dsW(p T - 8,q T ;p T ,q T ;T;i%(3)w(s,q T ) - 
^ J dr J dsW(p T ,q T ;p T - s,q T ;T;ihP)to(s,q T ), 



(8) 



where equation has to be supplemented by an initial condition for the spectral density, W(pi; p2', ih/3) = 
W(p,i; ii2 \ t — 0;ih/3), which can be expressed in the form of a finite difference approximation of the path inte- 
gral pi M, M ■ 
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with * {ni;H2;qi, ■■■,q n ;q' 1 , ...,q' n ;ihf3) 
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where (x\y) denotes the scalar product of two vectors x • y. In this expression the original (unknown) density matrix 
of the correlated system e~P( K+u ' has been decomposed into 2n factors, each at a 2n times higher temperature, 
with the inverse e = (3/2n and the corresponding high temperature DeBroglie wave length squared A 2 = 2irfi 2 e/m. 
This leads to a product of known high-temperature (weakly correlated) density matrices, however, at the price of 2n 
additional integrations over the intermediate coordinate vectors (over the "path"). This representation is exact in the 
limit n — > oo, and for finite n an error of order l/n occurs. 

Expression (^) has to be generalized to properly account for spin effects. This gives rise to an additional spin part of 
the density matrix, whereas exchange effects can be accounted for by the antisymmetrization of only one off-diagonal 

matrix element in (|l0|). As a result, in ( JlCj ) one matrix clement, for example (^q[ e~ e ^ q'^j has to be replaced by the 

antisymmctrized expression (q{ ' ' K 
single-particle density matrix 



q' 2 ) which can be written as the following sum of determinants involving the 

A 
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(12) 



Here Cfj = N\/s\(N — s)l, and \q[ a — q' 2 b \ denotes the distance between the vertices q[ and q' 2 of particles a and b, 
respectively. As a result of the summation over spin variables and all possible exchange permutations, the determinant 
carries a subscript s denoting the number of electrons having the same spin projection. To improve the accuracy of 
the obtained expression, in the total potential U being the sum of all pair interactions U a b, we will replace U a b — ► 
where U°f is the proper effective quantum pair potential, see below. For more details on the path integral concept, 
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we refer to Refs. ElL E2l E3l E4, pq], for recent applications of this approach to correlated Coulomb systems, cf. 



Let us now come back to the integral equation (g). One readily checks that its solution can be represented 
symbolically by an iteration series of the form 

W 1 = W t + K t T W T 

= W* + K l Ti W Tl + K^KgW* + K^K^K^W^ + ..., (13) 

where W l and W T1 are the initial quantum spectral densities evolving classically during time intervals [0, t] and [0, n], 
respectively, whereas K T l +1 are operators that govern the propagation from time r, to r<+i. 

Since the time correlation functions (Q) are linear functionals of the spectral density, for them the same series 
representation holds, 



c FA (t) = {2vh)- 2v J J d^dfn <Kmi; 1*2) W(jjh; /i 2 ;t; ihp) = (4>\w*) = 
(flw*) + {mt^) + {m^K\w^) + (^k^KIKIW^) + ... (14) 

where 0(/Ui; ^2) = F '(yUi)A(// 2 ) and the parentheses (. . . | . . .) denote integration over the phase spaces {/«i;/Lt 2 } as 
indicated in the first line of the equation. 

Note that the mean value F(t) of a quantum operator F can be represented in a form analogous to (p^): 

F(t) = (2nh)~ 2v f f dindfr IM+IM t; ihfi), (15) 



which allows us to increase the efficiency of the simulations. In the following Section we apply this scheme to a system 
of interacting electrons and random scatterers. 



III. QUANTUM DYNAMICS 

As an application, in this work we will consider a system composed of heavy particles (called scatterers) with mass 
m s and negatively charged electrons with mass m e . To avoid bound state effects due to attraction we consider in this 
case study only negatively charged scatterers, assuming a positve backgroud for charge neutrality. The influence of 
electron-scatterer attraction will be studied in a further publication. 

The possibility to convert a series like (Q into a form convenient for probabilistic interpretation allows us to 
apply Monte Carlo methods to its evaluation. According to the general theory of the Monte Carlo methods for 
solving linear integral equations, e.g. 130], one can simultaneously calculate all terms of the iteration series (LLj 



Using the basic ideas of 1 30 we have developed a Monte Carlo scheme, which provides domain sampling of the terms 



giving the main contribution to the series (|14|), cf. |L8|, 19 . This sampling also reduces the numerical expenses for 



calculating the integrals in each term. For simplicity, in this work, we take into account only the first term of iteration 



series (14), which is related to the propagation of the initial quantum distribution ( |ll| ) according to the Hamiltonian 
equation of motion. This term, however, does not describe pure classical dynamics but accounts for quantum effects 
|p if and, in fact, contains arbitrarily high powers of Planck's constant. The remaing terms of the iteration series 



describe momentum jumps [[L8| 19 32 which account for higher-order corrections to the classical dynamics of the 
quantum distribution ( pT[ ) which are expected to be relevant in the limit of high density. A detailed investigation 
of the conditions for which the contribution of the next terms of the iteration series should be taken into account is 
presented in ||^, |33j . 

The dynamical evolution of this system is studied along the trajectory pair in phase space, Eqs. (^|) and (^J), on time 
scales less than a maximum time t' , This time is chosen small enough such that the system of scatterers is practically 
stationary within the time interval < t < t' . The initial states of the system evolution are proper equilibrium states 
and thus relate to physical times much larger than t' . So the initial micro states (particle configurations in phase 
space) can be randomly generated and an ensemble averaging should be performed with the canonical density operator 
of the electron-scatterer system, where the latter has been computed by the path integral Monte-Carlo method NEW 
pjj according to the probability distribution, which is proportional to \^>\ (|l0|). So due to the time reversibility it 
is more convenient to start generation of the trajectories (Eqs. (|5|) and (^)) from time r = starting from sampled 
by Monte Carlo method initial particle configurations (points tj on FigT). The phase cofactor of the complex-valued 
function 'J was taken into account by introducing the weight function ]3^, |34], |35| of the initial configurations for the 
subsequent dynamic evolution. In principle, the method is also applicable to liquids or plasmas by assigning smaller 
masses to the scatterers. 
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This approach allows us to generate, in a controlled way, various kinds of initial conditions of the many-body system, 
in particular (i) those which are characteristic of the fully interacting system [i.e. including scatterer-scatterer (s-s), 
electron-scatterer (e-s), and electron-electron (e-e)] and (ii) those which result if some aspects of these interactions are 
ignored. In all cases, the short-time dynamics can then be followed by including or excluding the (e-e)-interaction. 

For the numerical calculations we introduce dimensionless units, using the maximum value of the time characteristic 
for the short-time dynamics, t' as the unit of the dynamic time t of the system. So the dimensionless time defined by 
6 = t/t 1 will always vary between and 1. As a unit of length we take the reciprocal wavenumber fc _1 , determined 
by the ratio k 2 — 2m e Eo/h 2 . E.g., for one electron in an external potential field U = VoU(q) the operator exponent 
of the time evolution propagator exp(—iHt/h) can be rewritten in the form 

where A is the Laplace operator, Vb is characteristic strength of the interaction potential, M = h/2Eot\ £o = 
V /2E M = Vot'/h. 

The tensor of the electrical conductivity is given by: 

a aj (uj)= J dte iuJt - et d\ (j^J a (t + ifi\)} (17) 
according to the quantum Kubo formula fl6|| , where e — * +0. The current operator is 

N e N e 



i=l 



j a =j2eqm=J2^ ( 18 ) 

III,- 



with q a being the a— component of the velocity operator of the electron. In the Wigner representation this tensor 
reads 



a ay (io)= ^ dte 1 ^ 1 f dX4> ai {t,X)- k2(7a ^ 
Jo Jo 



(19) 



e 2 2E Q t' 
where 

(j) ai (t,X) = ^J 7 J Q (i + ifrX)^ = (2nhy 2 ^ J d^id^ 2 J^(^i)J a (^2)W(^i;fi 2 ;t-,ih(3-,ihX), (20) 

and the spectral density W(/ii; H2',t;ihf3;ihX) is defined as in Eq. (|j|) by replacing t* by n = t% + ihX, and t c by 
T 2 = ^2 ~~ ih([3 — A). The dimensionless conductivity tensor is denoted by a a ^(u}). 

Our model of correlated electrons interacting with an array of random scatterers is described by the following 
Hamiltonian, 

H = H e + H es + H ss = (K e + U ee ) + U es + (K s + U„), 



11, = VAaVS'if , 

i=i 



\li - Qj\ 



N c N s 

i=l i=i 



Ae 



2Mm s ^ •' ^ V. A 

In the problem at hand we choose m e /m s ~ 1/2000 <C 1, therefore, on the time scale t' the scatterers are practically 
fixed. H ss is the Hamiltonian of the scatterers, and £q e = V Q ee t' /h, £g s = V es t' /Ti and £q s = V ss t' /Ti are the constants 
of the binary interaction between electrons (e-e) u (\qi — <7j|/A ee ), electrons and scatterers (e-s) u {\qi — Qj\/X es ) and 
between scatterers (s-s) u (\Qi — Qj\/X ss ), respectively, qi represent the positions of the electrons (i = 1, . . . , N e ), and 
Qj those of the scatterers (j = 1, . . . , N 3 ). 
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The (e-e), (e-s) and (s-s) interactions are taken in the form of the Kelbg potential |36|, |37j], which is the exact 
quantum pair potential of an ensemble of weakly nonidcal and weakly degenerate charged particles: 

1 r 2 1 

u(\x ab \) = 1 - e- % °» + ^TtXab (1 - erf (x ab )) (22) 

where x ab — \r a b\/^ab, with A being the thermal wave length given by \ 2 ab = Hab the reduced mass [i^ = 

m^ 1 + m^ 1 , V^ b = and e a , eb being the respective charges. The erf(x) stands for the error function. Note that 
the Kelbg potential is finite at zero distances and coincides with the Coulomb potential at distances larger than A a ^. 
This potential has recently been successfully applied to the computation of thermodynamic |25|, |2(| and dynamic 
properties |3£| of dense quantum plasmas by means of path integral Monte Carlo and classical molecular dynamics 
methods, respectively. It is, therefore, expected to provide an accurate description of the full Coulomb interaction 
between all particles within the quantum dynamics approach of the present paper. 

Due to the large difference in the masses of electrons and heavy scatterers we can use two simplifications. The 
first one consists in the antisymmetrization of the density matrix only for electron spins and space coordinates. The 
second one, as we mentioned above, is the use of the adiabatic approximation for the dynamical evolution, where the 
positions of the heavy scatterers in each initial configuration were fixed during the electron-dynamics time given by 
the scale t'. The dynamic evolution has been realized according to relations (|B|), (||). Let us stress, that the Kelbg 
potential appearing in the canonical density operator used for ensemble averaging in ( |To| ) has to be taken at the 
inverse temperature e = (3/2n while the simulation of the dynamic evolution, according to Eqs. (||||) involves the 
Kelbg potential at the temperature 1/(3. 



To simplify the computations we included in (12) only the dominant contribution to the sum over the total electron 
spin s corresponding to s = N/2 electrons having spin up and down, respectively. The contribution of the other terms 
decreases rapidly with particle number and vanishes in the thermodynamic limit. To further speed up the convergence 
of the numerical calculations we bounded the integration over the variables qM , Qi , Qm > Qi responsible for oscillations 
of the function ip in ( |TT| ) and checked the insensibility of the obtained results to these operations. 



IV. NUMERICAL RESULTS 



We now apply the numerical approach explained above to the problem of an interacting ensemble of electrons and 
disordered scatterers in one dimension. In all calculations and figures times, frequencies and distances are measured 
in atomic units. The average distance between electrons, R s = l/n e ao, was taken to be 12.0; 2.6; 0.55, the densities 
of electrons and heavy scatterers are equal. The results obtained were practically insensitive to the variation of the 
whole number of the particles in Monte Carlo cell from 30 upto 50 and also of the number of high temperature 
density matrices (determined by the number of factors n) in (|10|), ranging from 10 to 20. Estimates of the average 
statistical error gave the value of the order 5-7 percent. We studied two different temperatures: ksT /\Vq s \ = 0.45 
and 0.28, corresponding to X ee /ao ~ 2.2 and X ee /ao ~ 3.5, respectively. The strengths of the three interactions in the 
system are fixed arbitrarily at the ratio Vq c : Vq s : Vq s ~ 0.7 : 1 : 32. 

According to the Kubo formula jli| our calculations include two different stages: (i) the generation of the initial 
conditions (configuration of scatterers and electrons) in the canonical ensemble with probability proportional to the 
quantum density matrix and (ii) the generation of the dynamic trajectories on time scales t' in phase space starting 
from these initial configurations. 

The results presented below are related to three different situations. The first two situations refer to the generation 
of initial configurations used for ensemble averaging, where the (e-e) interaction is fully included (called "interacting 
ensemble"), while the electron dynamics was simulated in case I. with (e-e) interaction ("interacting dynamics") and, 
in case II., without (e-e) interactions ( "noninteracting dynamics"). For reference a third situation, III., was studied 
where the (e-e) interaction is neglected completely, i.e. both, in the initial ensembles ("noninteracting ensemble") 
and in the dynamical evolution, see also Fig. 1 for illustration. As we mentioned before to avoid the influence of (e-s) 
bound state effects we have considered only a system of negatively charged heavy scatterers here. The influence of 
(e-s) attraction on localization will be studied in a further publication. 



A. Temporal quantum momentum-momentum correlation functions 

First we discuss the influence of (e-e) interactions on the temporal quantum momentum-momentum correlation 
function. Fig. 2 shows momentum-momentum correlation functions corresponding to the situations I. and II. in- 
troduced above, i.e. to interacting ensembles with interacting (curves 1) and noninteracting (curves 2) dynamics, 
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respectively. Fig. 2 allows us to compare the obtained functions for two temperatures and densities varying over one 
order of magnitude. 

The traces clearly show coherent oscillations, which are a manifestation of localization tendencies even at finite 
temperature. Note that at higher temperatures (right-hand panel) these oscillations are less pronounced compared 
to lower temperatures (left-hand panel). The damping times of these coherent oscillations are clearly different for 
the noninteracting and interacting dynamics. We observe that the electron-electron interaction leads to a strong 
reduction (curve 1) of this damping time as compared to (e-e) noninteracting electrons (curve 2). Furthermore, for the 
noninteracting dynamics the appearance of deep aperiodic modulations is seen. The dynamics of (e-e) noninteracting 
electrons leads to at least one large oscillation and several damped small ones, reflecting a strong spatial confinement 
of the electron system in this case. For interacting dynamics both, the first aperiodic modulation and the subsequent 
damped coherent oscillations are less pronounced. 

The damping time for coherent oscillations increases with increasing density and decreasing temperature. The 
physical reason of this phenomenon is the tendency towards ordering of the scatterers, which is obvious from the 
corresponding pair distribution functions (g ss ) presented in Fig. 3, which develop periodic modulations for higher 
densities. 

As an interesting side-result Fig. 4 presents the influence of the choice of initial conditions. In the left column 
we show the momentum- momentum correlation functions corresponding to situations II. and III., i.e. interacting 
(curves 1) and noninteracting (curves 2) ensembles and, in both cases, dynamics without electron-electron interaction. 
The analysis of the equilibrium configurations and pair distribution functions for the noninteracting ensemble (right 
column of Fig. 4) shows that the main contributions to the ensemble averaging result from absolutely different 
particle configurations if compared to the case of the interacting ensembles discussed above. The (e-e) noninteracting 
electrons are gathered mainly in the extended deepest potential well, while the heavy scatterers form a potential 
barrier confining the electrons. The typical initial configuration for equilibrium averaging, therefore, contains a so 
called 'electron bubble'. The (e-c) noninteracting electrons oscillate inside this 'bubble' with frequencies defined by 
the eigenvalues of the resulting potential energy profile of this well. On the other hand, for the typical configurations 
of interacting ensembles (cf. Fig. 3) the scatterers and electrons are both distributed more or less uniformly in space 
and the potential profile shows only spatial small-scale random inhomogeneitics. 

In the case of noninteracting ensembles we can see high frequency coherent oscillations with frequencies defined by 
the potential forming the 'bubble', while for interacting ensembles with their spatial small-scale potential profile only 
aperiodic oscillations (cf. Fig. 2, curves 2) are seen. 

B. Fourier transform of the momentum— momentum time correlation functions 

Figures 5-7 present the real and imaginary parts of the diagonal elements of the electrical conductivity tensor 
versus frequency, i.e. the real and imaginary parts of the Fourier transform of the temporal momentum-momentum 
correlation functions. The real part of the Fourier transform characterizes the Ohmic absorption of electromagnetic 
energy and has the physical meaning of electron conductivity, while the imaginary part presents (e — 1)c<j, where e is 
the permittivity of the system. 

Curves 1 and 2 of Fig. 5 show an opposite behavior for vanishing frequency. Most remarkably, the low-frequency 
conductivity related to the interacting dynamics (curve 1) is positive, while that for the noninteracting dynamics 
(curve 2) has a maximum at some finite frequency related to the coherent oscillations in the time domain. For lower 
frequency it changes sign. However, the real part of the conductivity has to be nonzero. The negative contributions 
are due to weakly damped very slow oscillations with time scales exceeding the scale t' considered for the calculation of 
the dynamics for the noninteracting dynamics (see the temporal momentum-momentum correlation functions, Fig. 2). 
To overcome this deficiency of our model one has to increase the time t' and/or to take into account the slow motion 
of the heavy particles, which will destroy the coherent oscillation of the light electrons trapped by the heavy particles 
and thus suppress these negative values. In fact, we found that calculations performed for longer and longer times 
lead to decreasing negative contributions for low frequencies. 

Nevertheless, let us stress that the results obtained (Figs. 2 and 5) allow us to conclude that the behavior of the 
conductivity in the vicinity of zero frequency resembles the characteristic features of Anderson localization (i.e. a 
vanishing zero-frequency conductivity and enhanced oscillations related to the maximum in the real part at finite 
frequency) for noninteracting dynamics of electrons. The (e-e) interaction, on the other hand, leads to a strong 
increase of the conductivity at low frequencies and less pronounced oscillations. 

Let us note that the high frequency tails of the real and imaginary parts of the Fourier transforms presented in 
Figs. 5 and 6 coincide with each other, as the main contribution is due to the fast trajectories with high virtual energy 
This means, the (e-e) interaction does not influence the behavior of the high frequency tails. Furthermore, we have 
checked that the high- frequency tails of the real and imaginary parts of the Fourier transform can be described by 
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the Drude-asymptotics for free electrons. This reflects the fact that also the (e-s) interaction (its coupling constant 
being comparable to that of the (e-e) interaction) does not effect the behavior of the high-energy trajectories. 

The Fourier transform of the correlation functions for noninteracting dynamics and both interacting and nonin- 
teracting ensembles are presented in Fig. 7 (cf. also Fig. 4). The most remarkable feature of these figures is the 
splitting of the peak related to the coherent electron oscillations for the noninteracting ensembles. The positions of 
the inserted triangles show the energy difference between subsequent eigenenergy levels calculated for an averaged 
bubble. The comparison of the peak and the triangle positions confirms the fact that p6[ the Fourier transform of 
the momentum-momentum correlation function contains a sum of products of delta- functions and momentum matrix 
elements related to transitions between these energy levels. The imaginary part of the Fourier transform presented 
in the right column of Fig. 7 reflects that these characteristic features are related to the resonance oscillations in the 
'bubble' mentioned above. 

Figure 5 shows that for noninteracting dynamics, at low temperatures and moderate densities, e.g. R s — 2.6 (cf. 
Figs, b), a well-resolved absorption peak appears which is related to the electron energy level separation. Note that 
this peak is present even in the case of interacting ensembles. However, there, the height of this peak is considerably 
smaller than in the case of noninteracting ensembles (Fig. 7). 

C. Position dispersion 

Finally, we discuss the position dispersion of electrons for interacting and noninteracting ensembles and noninter- 
acting dynamics, cf. Fig. 8. Electron localization in the bubble results in a very slow growth (curve 2) of the position 
dispersion in comparison with the case of the small-scale potential profile (curve 1) of the initial configurations, as seen 
in Fig. 8. The square root of the position dispersion at times of the order of unity yields the typical size of the bubble. 
However, for larger times, t > 80, curve 2 shows a behavior, which is typical for particle diffusion. This characteristic 
change in behavior is due to the tail of the electron energy distribution function and is related to the fact that we are 
dealing with non-selfaveraging quantities. For the high-frequency tails of the Fourier transforms the virtual energy of 
the trajectories may be large. These exponentially rare fast trajectories can give an exponentially large contribution 
to the position dispersion, as the difference in positions of fast trajectories may be exponentially large at sufficiently 
large time. A similar problem connected with exponentially large contributions of exponentially rare events arises in 
the consideration of classical wave propagation in random media. There, it is well known ]3£| that in one-dimensional 
systems the dispersion of wave intensities is not self-averaging, as exponentially rare configurations of scatterers can 
give rise to exponentially large contributions of intensity at large distances from the wave source. 

V. CONCLUSION 

The Quantum-Dynamics-Monte-Carlo approach applied to a one-dimensional system of interacting electrons in an 
array of fixed random scatterers at finite temperature gives strong evidence for an enhancement of the quantum mobil- 
ity of electrons due to their mutual long-range many-particle interaction and thus substantiates previous expectations 
drawn from schematic models. 

As in our approach the temperature is taken to be nonzero, Anderson localization is not expected to show up in the 
results for noninteracting electrons in a strict sense. On the other hand, for the considered temperatures (comparable 
to or less than the coupling constant of the electron-scatterer potential) localization tendencies are clearly observable 
in this case. They manifest themselves both, in the low-frequency behavior of the momentum-momentum correlation 
function related to the conductivity of the electron system, and in coherent oscillations. For nonzero electron-electron 
interaction these localization tendencies are relaxed and the frequency- dependent correlation function has a form not 
unlike a Drude behavior. The high-frequency tails resemble a Drude-behavior of free electrons in all cases. 

For initial conditions describing a noninteracting ensemble, the scattering centers form electron-bubbles which 
accomodate the [(e-e) noninteracting] electrons in their eigenstates and lead to well-defined coherent oscillations. In 
contrast, for configurations describing interacting ensembles, the scatterers and electrons are both distributed more 
or less uniformly in space and the potential profile shows only spatial small-scale random inhomogeneities. Then 
the dynamics of interacting electrons in this spatial small-scale random potential profile shows only fast damped 
oscillations with much smaller damping time if compared to noninteracting dynamics and results in an increased 
low-frequency conductivity. 

The present work shows that it is possible to treat a system of many quantum particles interacting with each 
other via the long-range Coulomb potential in a numerically rigorous scheme. Here we have considered the heavy 
particles to be essentially immobile species. On the other hand, if their mass is being reduced they will take part 
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in the dynamics, and we then have a system resembling a plasma of, e.g., electrons and holes. Such calculations are 
presently under way and will be published shortly. 
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FIGURES 



Fig.l Illustration of the Quantum-Dynamics scheme for computation of the correlation function Cfa{t). The 
quantum general dynamics scheme generates particle trajectories (long arrows) in phase space p — q starting from 
an initial state at t = 0. The two arrows indicate different ensembles used for averaging over the initial states 
which do (do not) include e-e interaction: interacting ensemble, (noninteracting ensemble). At randomly chosen 
times t\,. . . ,ti, t i+ i, . . . ,tx, a pair of dynamic trajectories (short arrows) is propagated in positive and negative time 
direction, e.g. from ti to ti + t/2 and t{ — r/2, respectively. Cfa(t) is averaged over all K trajectories. Bold (dashed) 
arrows indicate (non-)interacting dynamics including (neglecting) e-e interaction. 

Fig. 2 Temporal momentum-momentum correlation for interacting ensembles for the case of interacting (1) and 
noninteracting (2) dynamics. Temperatures are I<;bT/\Vq s \ = 0.28 (left panel: a,b) and k,BT/\VQ S \ — 0.45 (right 
panel: c,d). 

Fig. 3 e-e, e-s and s-s pair distribution functions for the interacting ensemble, for the same parameters as in Fig. 2. 

Fig. 4 Left Column: Temporal momentum-momentum correlation functions for the case of noninteracting dynamics 
and interacting (1) and noninteracting (2) ensembles. Temperature is fcsT/|V^f s | = 0.45. Right Column: Pair 
distribution functions for the noninteracting ensembles. 

Fig. 5 Real part of the Fourier transform of the temporal momentum-momentum correlation functions of Fig. 2 

Fig. 6 Imaginary part of the Fourier transform of the temporal momentum-momentum correlation functions of Fig. 

2. 

Fig. 7 Real part (left column) and imaginary part (right column) of the Fourier transform of the temporal momentum- 
momentum correlation functions of Fig. 4. Triangles indicate the transition energies between the lowest energy 
eigenvalues for an anlytical model for the deepest scatterer potential which traps the electrons (vertical coordinates 
of the triangles are arbitrary). 

Fig. 8 Electron position dispersion for noninteracting dynamics and interacting (1) versus noninteracting (2) ensem- 
bles, cf. Fig. 4. 
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